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^ : ABSTRACT 

As an explosion develops in the collapsed core of a massive star, neutrino emission drives 
convection in a hot bubble of radiation, nucleons, and pairs just outside a proto-neutron star. 
Shortly thereafter, neutrinos drive a wind-like outflow from the neutron star. In both the convec- 
tive bubble and the early wind, weak interactions temporarily cause a proton excess (Y c &0.50) 
to develop in the ejected matter. This situation lasts for at least the first second, and the ap- 
proximately 0.05 - 0.1 Mq that is ejected has an unusual composition that may be important for 
nucleosynthesis. Using tracer particles to follow the conditions in a two-dimensional model of a 
successful supernova explosion calculated by Janka, Buras, & Rampp (2003), we determine the 
composition of this material. Most of it is helium and 56 Ni. The rest is relatively rare species 
produced by the decay of proton-rich isotopes unstable to positron emission. In the absence of 
pronounced charged-current neutrino capture, nuclear flow will be held up by long-lived waiting 
point nuclei in the vicinity of 64 Ge. The resulting abundance pattern can be modestly rich in a 
few interesting rare isotopes like 45 Sc, 49 Ti, and 64 Zn. The present calculations imply yields that, 
when compared with the production of major species in the rest of the supernova, are about those 
needed to account for the solar abundance of 45 Sc and 49 Ti. Since the synthesis will be nearly 
the same in stars of high and low metallicity, the primary production of these species may have 
discernible signatures in the abundances of low metallicity stars. We also discuss uncertainties 
in the nuclear physics and early supernova evolution to which abundances of interesting nuclei 
are sensitive. 



Subject headings: supernovae, nucleosynthesis 



1. INTRODUCTION 

When the iron core of a massive star collapses 
to a neutron star, a hot proto-neutron star is 
formed which radiates away its final binding en- 
ergy as neutrinos. Interaction of these neutrinos 
with the infalling matter has long been thought 
to be the mechanism responsible for exploding 
that part of the progenitor external to the neu- 
tron star and making a supernova (e.g., Janka 
2001; Woosley, Heger, & Weaver 2002, and ref- 
erences therein). During the few tenths of a sec- 
ond when the explosion is developing, a convec- 
tive bubble of photo-disintegrated matter (nucle- 
ons), radiation, and pairs lies above the neutron 
star but beneath an accretion shock. Neutrino 
interactions in this bubble power its expansion, 
drive convective overturn, and determine its com- 
position. Since baryons exist in the bubble only 
as nucleons, the critical quantity for nucleosyn- 
thesis is the proton mass fraction (Y c ). Initially, 
in part because of an excess of electron neutri- 
nos over antineutrinos, Y e z0.5 (Qian & Woosley 
1996). As time passes, however, the fluxes of the 
different neutrino flavors and their spectra change 
so that Y c evolves and becomes considerably less 
than 0.5. This epoch, also known as the "neutrino- 
powered wind", has been explored extensively as 
a possible site for the r-process (Qian & Woosley 
1996; Hoffman, Woosley, & Qian 1997; Woosley et 
al. 1994; Cardall & Fuller 1997; Qian & Wasser- 
burg 2000; Takahashi, Witti, & Janka 1994; Ot- 
suki et al. 2000; Sumiyoshi et al. 2000; Thompson, 
Burrows, & Meyer 2001) as well as 64 Zn and some 
light p-process nuclei (Hoffman et al. 1996). 

In this paper we consider nucleosynthesis dur- 
ing the earlier epoch when Y e is still greater than 
0.5. This results in a novel situation in which the 
alpha-rich freeze out occurs in the presence of a 
non-trivial abundance of free protons. The re- 
sulting nuclear flows thus have characteristics of 
both the alpha-rich freeze out (Woosley, Arnett, 
& Clayton 1973; Woosley & Hoffman 1992) and 
the rp-process (Wallace & Woosley 1981). Several 
proton-rich nuclei, e.g., 64 Ge and 45 Cr, are pro- 
duced in such great abundance that, after ejection 
and decay, they contribute a significant fraction of 
the solar inventory of such species. 



2. Supernova Model and Nuclear Physics 
Employed 

2.1. Explosion Model for a 15 M Star 

The nucleosynthesis calculations in this paper 
are based on a simulation of the neutrino-driven 
explosion of a nonrotating 15Mq star (Model 
S15A of Woosley & Weaver 1995) by Janka, Buras, 
& Rampp (2003) (see also Janka et al. 2004). The 
post-bounce evolution of the model was followed 
in two dimensions (2D) with a polar coordinate 
grid of 400 (noncquidistant) radial zones and 32 
lateral zones (2.7 degrees resolution), assuming az- 
imuthal symmetry and using periodic conditions 
at the boundaries of the lateral wedge at ±43.2° 
above and below the equatorial plane. Convection 
was seeded in this simulation by velocity pertur- 
bations of order 10 -3 , imposed randomly on the 
spherical post-bounce core. 

The neutrino transport was decribed by solv- 
ing the energy-dependent neutrino number, en- 
ergy, and momentum equations in radial direction 
in all angular bins of the grid, using closure rela- 
tions from a model Boltzmann equation (Rampp 
& Janka 2002). Neutrino pressure gradients and 
neutrino advection in lateral direction were taken 
into account (for details, see Buras et al. 2004). 
General relativistic effects were approximately in- 
cluded as described by Rampp & Janka (2002). 

Although convective activity develops in the 
neutrino-heating layer behind the supernova (SN) 
shock on a time scale of several ten milliseconds 
after bounce, no explosions were obtained with the 
described setup until ~250ms (Buras et al. 2003), 
at which time the very CPU-intense simulations 
usually had to be terminated. The explosion in 
the simulation discussed here was a consequence 
of omitting the velocity-dependent terms from the 
neutrino momentum equation. This manipulation 
increased the neutrino-energy density und thus the 
neutrino energy deposition in the heating region 
by ^20-30% and was sufficient to convert a failed 
model into an exploding one (see also Janka et al. 
2004, Buras et al. 2004). This sensitivity of the 
outcome of the simulation to only modest changes 
of the transport treatment demonstrates how close 
the convecting, 2D models of Buras et al. (2003) 
with energy-dependent neutrino transport are to 
ultimate success. 
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The evolution from the onset of core collapse 
(at about —175 ms) through core bounce and con- 
vective phase to explosion is shown in terms of 
mass shell trajectories in Fig. 1. The explosion 
sets in when the infalling interface between Si layer 
and oxygen-enriched Si layer reaches the shock 
at about 160 ms post bounce. The corresponding 
steep drop of the density and mass accretion rate, 
associated with an entropy increase by a factor of 
~ 2, allow the shock to expand and convection to 
become more violent, thus establishing runaway 
conditions. The calculation was performed in 2D 
for following the ejection of the convective shell 
until 470 ms after bounce. While matter is chan- 
neled in narrow downflows towards the gain ra- 
dius, where it is heated by neutrinos and some of 
it starts expanding again in high-entropy bubbles, 
its neutron-to-proton ratio is set by weak interac- 
tions with electron neutrinos and antineutrinos as 
well as electron and positron captures on free nu- 
cleons. The final value of Y e is a crucial parameter 
for the subsequent nucleosynthesis. The mass dis- 
tribution of neutrino-heated and -processed ejecta 
from the convective bubble is plotted in Fig. 2. 

At 470 ms after bounce the model was mapped 
to a ID grid and the subsequent evolution was sim- 
ulated until 1300 ms after bounce. With accretion 
flows to the neutron star having ceased, this phase 
is characterized by an essentially spherically sym- 
metric outflow of matter from the nascent neutron 
star, which is driven by neutrino-energy deposi- 
tion outside the neutrinosphere (Woosley & Baron 
1992; Duncan, Shapiro, & Wasserman 1986). 
This neutrino-powered wind is visible in Fig. 1 
after ^500 ms. The fast wind collides with the 
dense shell of slower ejecta behind the shock and 
is decelerated again. The corresponding negative 
velocity gradient steepens to a reverse shock when 
the wind expansion becomes supersonic (Fig. 1; 
Janka & Muller 1995). Characteristic parameters 
for some mass shells in this early wind phase are 
shown in Fig. 3. Six representative shells are suf- 
ficient, because the differences between the shells 
evolve slowly with time according to the slow vari- 
ation of the conditions (neutron star radius, gravi- 
tational potential, neutrino luminosities and spec- 
tra) in the driving region of the wind near the neu- 
tron star surface. In Table 3 the masses associated 
with the different shells are listed. 

At the end of the simulated evolution the model 



has accumulated an explosion energy of approxi- 
mately 0.6 x 10 51 erg. The mass cut and thus initial 
baryonic mass of the neutron star is 1.41 M Q . The 
model fulfills fundamental constraints for Type II 
SN nucleosynthesis (Hoffman et al. 1996) because 
the ejected mass having Y e < 0.47 is <10~ 4 M Q 
(see Fig. 2) and thus the overproduction of N=50 
(closed neutron shell) nuclei of previous explosion 
models does not occur. More than 83% of the 
ejected mass in the convective bubble and early 
wind phase (in total 0.03 M Q in this rather low- 
energetic explosion) have Y c > 0.5. The ejec- 
tion of mostly p-rich matter is in agreement with 
ID general relativistic SN simulations with Boltz- 
mann neutrino transport in which the explosion 
was launched by artificially enhancing the neu- 
trino energy deposition in the gain layer (Thiele- 
mann et al. 2003; Frohlich et al. 2004). The 
reason for the proton excess is the capture of elec- 
tron neutrinos and positrons on neutrons, which is 
favored relative to the inverse reactions because of 
the mass difference between neutrons and protons 
and because electron degeneracy becomes negligi- 
ble in the neutrino-heated ejecta (Frohlich et al. 
2004; Qian & Woosley 1996). 

Although the explosion in the considered SN 
model of Janka, Buras, & Rampp (2003) was 
obtained by a regression from the most accu- 
rate treatment of the neutrino transport, it not 
only demonstrates the proximity of such accurate 
models to explosions, but also provides a con- 
sistent description of the onset of the SN explo- 
sion due to the convectively supported neutrino- 
heating mechanism, and of the early SN evolution. 
The properties of the resulting explosion are very 
interesting, including the conditions for nucleosyn- 
thesis. The Y e values of the ejecta should be rather 
insensitive to the manipulation which enabled the 
explosion. On the one hand the expansion veloc- 
ities of the high-entropy ejecta are still fairly low 
(less than a few 10 8 cms _1 ) when weak interac- 
tions freeze out, and on the other hand the omit- 
ted velocity-dependent effects affect neutrinos and 
antineutrinos in the same way. 

2.2. Outflows in the Convective Bubble 

In order to calculate the nucleosynthesis it is 
necessary to have a starting composition and the 
temperature-density (T — p) history of the matter 
as it expands and is ejected from the supernova. 
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Because the matter is initially in nuclear statistical 
equilibrium, the initial values of Y e , T, and p deter- 
mine the composition which is just protons with a 
mass fraction Y e and neutrons. We are most inter- 
ested in the innermost few hundredths to one tenth 
of a solar mass to be ejected. This matter has an 
interesting history. It was initially part of the sil- 
icon shell of the star, but fell in when the core 
collapsed, passed through the SN shock and was 
photodisintcgrated to nucleons. Neutrino heating 
then raised the entropy and energy of the matter 
causing it to convect. Eventually some portion of 
this matter gained enough energy to expand and 
escape from the neutron star, pushing ahead of it 
the rest of the star. As it cooled, the nucleons 
reassembled first into helium and then into heavy 
elements. 

The temperature-density history of such matter 
is thus not given by the simple ansatz often em- 
ployed in explosive nucleosynthesis — "adiabatic 
expansion on a hydrodynamic time scale" . In fact, 
owing to convection, the temperature history may 
not even be monotonic. Here we rely on tracer 
particles embedded in the so called "hot connec- 
tive bubble" of the 15 M Q SN model calculated by 
Janka, Buras, & Rampp (2003) (Fig. 4). These 
tracer particles were not distributed uniformly in 
mass, but chosen to represent a range of Y e in the 
ejecta. 

The proton-rich outflows of interest here be- 
gin at about 190 ms after core bounce (Fig. 1). 
Entropies and electron fractions characteristic of 
a few different trajectories are given in Table 2. 
Each trajectory represents a different mass ele- 
ment in the convective bubble. As is seen, Y c 
for the different trajectories lies in the range from 
0.5 — 0.546, and the entropies per nucleon are mod- 
est, s/ki, ~ 30 — 50. Figure 2 shows the ejected 
mass versus Y e during the convective phase of the 
SN explosion. 

At the end of the 2D calculation of Janka, 
Buras, & Rampp (2003), the mass element in a 
typical trajectory had reached a radius of about 
2000 km (corresponding to the time when the 
SN model was mapped from 2D to ID and thus 
detailed information for the mass elements was 
lost). Temperatures at this radius were typically 
T g = T/10 9 K « 4-5, which is still hot enough 
that nuclei have not yet completely re-assembled. 
To follow the nucleosynthesis until all nuclear re- 



actions had frozen out it was necessary to extrap- 
olate the trajectories to low temperature. In doing 
so, we assumed that the electron fraction and en- 
tropy were constant during the extrapolated por- 
tion of the trajectory. This should be valid because 
the number of neutrino captures suffered by nuclei 
beyond ~ 2000 km is small. 

We considered two approximations to the ex- 
pansion which should bracket the actual behav- 
ior. The first assumes homologous expansion at 
a velocity given by the Janka et al. calculation 
between 10 billion and 4 billion K. This ignores 
any deceleration experienced as the hot bubble en- 
counters the overlying star and is surely an under- 
estimate of the actual cooling time (though per- 
haps realistic for the accretion-induced collapse of 
a bare white dwarf). In particular, we estimated 
the homologous expansion time scale for each tra- 
jectory as Thom = (if — ^i)/l n (Pi/Vf) where the 
subscript i denotes the value of a quantity when 
T 9 = 10 and the subscript f denotes the value of a 
quantity at the last time given for the tracer par- 
ticle history (if « 436 ms, Tgj « 4-5). Values of 
Thom for different trajectories are given in Table 2. 

The second approximation was an attempt to 
realistically represent material catching up with 
the supernova shock. This extrapolation is based 
on smoothly merging the trajectories found in the 
calculations of Janka et al. with those calculated 
for the inner zone of the same 15 Mq supernova 
by Woosley & Weaver (1995). There are some 
differences. The earlier study was in one dimen- 
sion and the shock was launched artificially us- 
ing a piston. The kinetic energy at infinity of the 
Woosley- Weaver model was 1.2 x 10 51 erg; that of 
the Janka et al. model was 0.6 x 10 51 erg. Still 
the calculations agreed roughly in the tempera- 
ture and density at the time when the evaluation 
of tracer particles in the current 2D simulation was 
stopped. In order not to have discontinuities in the 
entropy at the time when the two calculations are 
matched, the density in the previous ID calcula- 
tion is changed slightly. This merging of the late 
time trajectories is expected to be reasonable be- 
cause the shock evolution at several seconds post 
core bounce is determined mostly by the explosion 
energy. 

We shall see in Sect. 3 that abundances of key 
nuclei are particularly sensitive to the time it takes 
the flow to cool from 2 • 10 9 K to 1 • 10 9 K. The ho- 
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mologous expansion approximation gives this time 
as about 100-200 ms, while the Kepler based es- 
timate gives this time as about 1 sec. Both esti- 
mates are rough and should be viewed as repre- 
senting upper and lowed bounds to the time scale. 

Figure 5 shows the evolution of density in a 
representative trajectory for each of the two ap- 
proximations to the flow at large radii. The tem- 
perature history in these trajectories is shown in 
Fig. 6. Note the irregular and non-monotonic evo- 
lution of the thermodynamic quantities at early 
times. 

2.3. Outflows in the Early Wind 

While the shock sweeps through and expels 
the stellar mantle, matter is still being contin- 
uously ablated from the surface of the cooling 
neutron star. Neutrino heating, principally via 
charged current neutrino capture, acts to maintain 
pressure-driven outflow in the tenuous atmosphere 
formed by the ablated material. This outflow has 
a higher entropy and is less irregular than the con- 
vective bubble. 

The evolution of material at radii smaller than 
a few hundred km is set by characteristics of the 
cooling neutron star, ft is at these small radii that 
the asymptotic entropy s and electron fraction Y e 
are set. At early times the neutron star has yet 
to radiate away the bulk of its gravitational en- 
ergy and so has a relatively large radius. Material 
escaping the star during this period only needs 
to gain a little energy through heating to escape 
the still shallow gravitational potential. Conse- 
quently, the entropy of the asymptotic outflow is 
about a factor of two smaller than the entropy of 
winds leaving the neutron star ~10 seconds post 
core-bounce. This can be seen from the analytic 
estimate provided by Qian & Woosley (1996) 

S « 235(L Pe , 51 4,Mev)- 1/6 ( 1 ^) 2/3 - (1) 

Here £p e , 5 i = £ Pe /10 51 erg/sec , e Pe ,McV approxi- 
mately the mean energy of electron anti-neutrinos 
and R is the neutron star radius. A lower entropy 
implies a higher density and therefore faster par- 
ticle capture rates at a given temperature. For 
proton-rich outflows this typically results in syn- 
thesis of heavier elements. 

The electron fraction in the outflow is set by a 



competition between different lepton capture pro- 
cesses on free nucleons: 

v e + n < > p + eT , (2) 

e + + n < > p + v e - (3) 

Because the neutron star is still deleptonizing at 
early times, the i/ e and v e spectra can be quite 
similar. Also, once heating raises the entropy of 
material leaving the neutron star, the number den- 
sities and spectra of electrons and positrons within 
the material become similar. Under these circum- 
stances the 1.29 MeV threshold for p — > n results 
in P e /e~ capture rates which are slower than the 
inverse v e /e + capture rates. Weak processes then 
drive the outflow proton rich. The electron frac- 
tion in the wind is mostly set by the competition 
between u e and D e capture (because captures 
freeze out when the density and temperature in 
the outflow become low, whereas high-energy neu- 
trinos streaming out from the neutrinosphere still 
continue to react with nucleons). When the com- 
position comes to equilibrium with the neutrino 
fluxes, 



Here \ v represents the electron neutrino or an- 
tineutrino capture rate on neutrons or protons. 
Because the star is still deleptonizing at early 
times, the v e and v e spectra can be quite similar. 
The 1.29 MeV threshold for v e capture then leads 
to A„ e „ > A PeP , and proton-richness is established 
in the outflow. Finally also the neutrino reactions 
cease because of the 1 /r 2 dilution of the neutrino 
density with growing distance r from the neutron 
star. 

Table 3 gives characteristics of the early wind 
found in the simulations of Janka, Buras, & 
Rampp (2003). As expected, the wind is pro- 
ton rich at early times. Eventually, the hardening 
of the v e spectrum relative to the v e spectrum 
will cause Y e to fall below 1/2. This turnover has 
not yet occurred when the hydrodynamic simula- 
tion was stopped. It should take place at a later 
time when the wind properties (mass loss rate, 
entropy) have changed such that the nucleosyn- 
thesis constraints for the amount of Y e < 0.47 
ejecta (Hoffman et al. 1996) will not be violated. 
At 1.3 s after bounce the mass loss rate of about 
3 x IO^MqS" 1 and wind entropy of ~ 80 fc b per 
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nucleoli in the Janka et al. model are likely to still 
cause an overproduction of N=50 nuclei if Y e went 
significantly below 0.5. 

The temperature in the wind at the end of the 
traced shell expansion is T 9 w 2 (Fig. 3). Approx- 
imations for the wind evolution at lower tempera- 
tures are the same as those discussed above. 

2.4. Nuclear Physics Employed 

The reaction network used for the present cal- 
culations is given in Table 1 . Estimates of reaction 
rates and nuclear properties used in our calcula- 
tions are the same as those used in the study of 
X-ray bursts by Woosley et al. (2004). Briefly, re- 
action rates were taken from experiment whenever 
possible, from detailed shell-model based calcula- 
tions (Fisker et al. 2001) for a few key (p, 7) rates, 
and from Hauser-Feshbach calculations (Rauscher 
& Thielemann 2000) otherwise. Proton separa- 
tion energies, which are crucial determinants of 
nucleosynthesis in flows with Y e > 1/2, were taken 
from a combination of experiment (Audi & Wap- 
stra 1995), the Hartree-Fock Coulomb displace- 
ment calculations of Brown et al. (2002) for many 
important nuclei with Z>N, and theoretical esti- 
mates (Moller et al. 1995). Choosing the best nu- 
clear binding energies is somewhat involved and 
we refer the reader to the discussion in Brown et 
al. (2002) and Fig. 1 of Woosley et al. (2004). 
Ground-state weak lifetimes are experimentally 
well determined for the nuclei important in this 
paper. At temperatures larger than 10 9 K the in- 
fluence of thermal effects on weak decays was es- 
timated from the compilation of Fuller, Fowler, & 
Newman (1982) where available. Table 7 gives the 
nuclei for which the Fuller et al. rates were used. 
A test calculation in which we switched thermal 
rates off and used only experimentally determined 
ground-state rates showed little effect on the im- 
portant abundances. Section 3.1 contains a dis- 
cussion of the influence of nuclear uncertainties on 
yields of some interesting nuclei. 

3. Nucleosynthesis Results 

Table 2 gives the major calculated production 
factors for a number of trajectories in the convec- 
tive bubble and for our two different estimates of 
the material expansion rate at low temperatures. 
Table 3 gives production factors for nuclei syn- 



thesized in different mass elements comprising the 
early wind. Here the production factor for nuclide 
i is defined as 

Pj = ^a|? (5) 

where M is the total mass in a given trajectory, 
M Ci — 13.5 M Q is the total mass ejected in the 
SN explosion, Xi is the mass fraction of nuclide 
i in the trajectory, and Xqj is the mass fraction 
of nuclide i in the sun. To aid in interpreting the 
tables we show in Fig. 7 plots of X^/Xq^ charac- 
terizing the nucleosynthesis in two representative 
hot-bubble trajectories. 

Production factors integrated over the differ- 
ent bubble trajectories are given in Table 4. If 
one assumes rapid expansion, production factors 
of 45 Sc, 63 Cu, 49 Ti, and 59 Co are all above 1.5. 
For the slower expansion time scale below 4 x 10 9 
K, which we regard as more realistic, a differ- 
ent set of nuclei are produced, especially 49 Ti and 
64 Zn. Depending upon mass and mctallicity, 49 Ti 
may already be well produced in other regions 
of the same supernova (Woosley & Weaver 1995; 
Rauscher, Hcger, Hoffman, & Woosley 2002), but 
64 Zn is not. The synthesis here thus represents a 
new way of making 64 Zn and this same process will 
function as well in zero and low metallicity stars as 
in supernovae today. However, 64 Zn was already 
known to be produced, probably in greater quan- 
tities, by the neutrino-powered wind (Hoffman et 
al. 1996). 

Production factors integrated over the different 
wind trajectories are given in Table 5. The some- 
what high-entropy wind synthesizes 45 Sc, 49 Ti and 
46 Ti more efficiently than the bubble. Typical val- 
ues of X/Xq for these three nuclei are approxi- 
mately 10 4 in the wind and approximately 2 • 10 3 
in the bubble. In the present calculations the in- 
tegrated production factor for Sc in the wind is 
between about 1.5 and 4.7 depending on the time 
scale describing the wind expansion at Tg < 2. 

For comparison, in the 15 M Q supernova of 
Rauscher, Heger, Hoffman, & Woosley (2002), this 
production factor was about 7 for many major 
species, including oxygen. This is close to the 
combined wind/bubble production factors of Sc 
and 46 > 49 Ti in the present calculations. The other 
most abundant productions in Tables 4 and 5 fall 
short of this - but not by much. The bulk pro- 
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duction factors in a 25 Mq supernova are about 
twice those in a 15, but our explosion model is not 
easily extrapolable to stars of other masses. If 25 
M Q stars explode with a similar kinetic energy it 
will probably take a more powerful central engine 
to overcome their greater binding energy and ac- 
cretion rate during the explosion. Probably this 
requires more mass in the convective bubble. In 
fact, the energy of the 15 M supernova used here, 
0.6 x 10 51 erg, would be regarded by many as low. 
It may be that the mass here should be doubled 
too. 

It is important to note that the species listed 
in Tables 4 and 5 are not made as themselves but 
as proton-rich radioactive progenitors. Major pro- 
genitors of important product nuclei are given in 
the far right column of Table 4. Typical progeni- 
tors of important nuclei are 3-4 charge units from 
stability. This can be understood through con- 
sideration of the Saha equation. Before charged 
particle reactions freeze out at T 9 w 1.5 — 2, nu- 
clear abundances along an isotonic chain are well 
approximated as being in local statistical equilib- 



rium: 



A(Z + 1,N) 
X(Z,N) 



in-5 /t\ P 5 ^Z+1>N 

10 exp(5 p /T) . 



'Z,N 



(6) 

Here S p is the proton separation energy of the 
Z+1,N nuclide, G represents the partition func- 
tion, p 5 = /Vl0 5 gcm- 3 , T 9 = T/10 9 K, and 
A=Z+N. Equation (6) predicts that the abun- 
dances of nuclei with S p < 500 keV are very small. 

Perhaps the most notable feature of the proton- 
rich trajectories is their inefficiency at synthesiz- 
ing elements with A>60. Neutron-rich outflows, 
by contrast, readily synthesize nuclides with mass 
A ~ 100. This is shown in Table 6 which gives 
production factors characterizing nucleosynthesis 
in somewhat neutron-rich winds occurring in the 
SN. The Kepler-based extrapolation of the first 
trajectory in Table 2 is used for these Y c < 0.5 
calculations. Estimates of the mass in each Y e bin 
for the calculations of Janka, Buras, & Rampp 
(2003) are shown in Fig. 2. 

Termination of the nuclear flow at low mass 
number in proton-rich outflows has a simple ex- 
planation. Unlike nuclei at the neutron drip 
lines, proton-rich waiting point nuclei have life- 
times much longer than the time scales charac- 



terizing expansion of neutrino-driven outflows. In 
addition, proton capture from waiting point nu- 
clei to more rapidly decaying nuclei is inefficient. 
To illustrate the difficulty with rapidly assembling 
heavier proton-rich nuclei, consider nuclear flow 
through 64 Ge. This waiting point nucleus has a 
lifetime of approximately 64 sec. The ratio of the 
amount of flow leaving 65 As to that leaving 64 Ge 
is found from application of the Saha equation 
above, 



A + ( 65 As)y( 65 As) 
A+( 64 Ge)y( 64 Ge) 



10- 2 -§^exp(S p /T). (7) 



Here A + represents the [3 + decay rate and S p is the 
proton separation energy of As. For 65 As, A + « 
ln(2)/0.1scc and for 64 Ge, A+ « ln(2)/64scc. By 
definition, proton capture daughters of waiting 
point nuclei are characterized by small proton sep- 
aration energies. The binding energy of 65 As still 
has large uncertainties, though is known to be less 
than about 200kcV (Brown ct al. 2002). Positron 
decay out of the proton capture daughter of the 
waiting point nuclei is negligible for such small 
proton separation energies. These considerations 
do not hold for X-ray bursts, where time scales 
characterizing nuclear burning can be tens or hun- 
dreds of seconds. 

The difficulty with rapid assembly of heavy 
proton-rich nuclei is also evident in the final free 
proton and alpha particle mass fractions. The 
trend of X p and X a with Y c is shown in Fig. 8 for 
the different Kepler extrapolated bubble trajecto- 
ries. Also shown in this figure is the proton mass 
fraction calculated under the assumption that all 
available nucleons are bound into alpha particles. 
This is an approximate measure of the mass frac- 
tion of available protons. Note that the mass frac- 
tion of protons in the two calculations are nearly 
identical. This is because assembly of proton-rich 
nuclei occurs on a very slow time scale set by a 
few f3 + rates. 

Because nucleosynthesis past A ~ 60 is ineffi- 
cient these proton-rich flows do not produce N=50 
closed shell nuclei. Historically, overproduction 
of N=50 nuclei has plagued calculations of super- 
nova nucleosynthesis (Howard et al. 1993; Witti et 
al. 1993; Woosley et al. 1994). The influence of 
weak interactions in driving some of the outflow 
to Y c > 1/2 ameliorates this problem. 
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3.1. Details of the Nucleosynthesis and 
Critical Nuclear Physics 

To aid in understanding the general character 
of these proton-rich flows we show in Fig. 9 the 
evolution of nuclear mass fractions as a function 
of the neutron number. At T 9 « 4, a captures 
have led to efficient synthesis of tightly bound 
species with N=28 and N=30. As temperature 
decreases a capture becomes less efficient and /3 + 
decay drives flow to higher neutron number. From 
Table 4 it is seen that the nuclei we are most inter- 
ested in arise from decay of nuclei with N=21, 24, 
31 and 32. From Fig. 9 it is clear that synthesis 
of nuclei with these neutron numbers represents 
a minor perturbation on the nucleosynthesis as a 
whole. 

Tables 4 and 5 show that 45 Sc, the only stable 
scandium isotope, has a combined wind/bubble 
production factor of about 6 if freeze-out is rapid 
and a combined production factor about 50% 
smaller in the slower Kepler extrapolated trajec- 
tories. Efficient synthesis of scandium in proton 
rich outflows associated with Gamma Ray Bursts 
has been noted previously by Pruet et al. (2004) , 
while Maeda & Nomoto (2003) found that scan- 
dium may also be synthesized explosively in shocks 
exploding anomalously energetic supernovac. In- 
deed, values presented here for Y e , s/kb, and r in 
the early SN wind are very close to estimates of 
these quantities in winds leaving the inner regions 
of accretion disks powering collapsars (MacFadyen 
& Woosley 1999; Pruet et al. 2004). The origin of 
Sc is currently uncertain and it may be quite abun- 
dant in low metallicity stars (Cayrel et al. 2004) 
suggesting a primary origin. In the present calcu- 
lations the yields of this element are close to those 
needed to explain the current inventory of Sc. 

To understand how synthesis of scandium de- 
pends on the outflow parameters and nuclear 
physics, note that Sc arises mostly from (3 + decay 
originating with the quasi waiting-point nucleus 
45 Cr. In turn, N=21 isotones of 45 Cr originate 
from (3 + decay out of isotones of 40 Ca. The dou- 
bly magic nucleus 40 Ca is efficiently synthesized 
through a sequence of alpha captures. At tem- 
peratures larger than about 2 • 10 9 K statistical 
equilibrium keeps almost all N=20 nuclei locked 
into 40 Ca. This nucleus is {3 stable and has a first 
excited state at 3.3 MeV, too high to be thermally 



populated. Flow out of N=20 can only proceed 
when the temperature drops to approximately 1.5 
billion degrees and statistical equilibrium favors 
population of 42 Ti over 40 Ca. The proton capture 
daughter of 40 Ca ( 41 Sc) has a proton separation 
energy of only 1.7 MeV and is not appreciably 
abundant. Decay out of 42 Ti is then responsible 
for allowing flow to N=21. 42 Ti has a well deter- 
mined j3 + half life of 199±6 ms, a proton sepa- 
ration energy which is uncertain only by about 5 
keV, and a first excited state too high in excitation 
energy to play a role in allowing flow to N=21. In 
short, nuclear properties are well determined for 
important N=20 nuclei. Once nuclei make their 
way to N=21 at Tg « 1.5, their abundances are 
divided between the tightly bound 45 Cr and 43 Ti. 
Here uncertainties in nuclear physics may be more 
important. For 45 Cr the proton separation energy 
is uncertain to about 100 keV and the spin of the 
ground state is uncertain. To the extent that the 
relative abundances are set by the Saha equation, 
these uncertainties could imply an uncertainty of 
a factor of several in the relative abundances of 
45 Cr and 45 Ti at T 9 « 1.5. In turn, this implies 
appreciable uncertainty in the estimated Sc yield. 

Whether or not Sc is efficiently synthesized fol- 
lowing decay of 45 Cr depends on the expansion 
time scale at low temperatures. This is because 
the (3 + daughter of 45 Cr is 45 V, which has a rela- 
tively small proton separation energy of 1.6 MeV. 
At low temperatures the Saha equation favors pro- 
ton capture to 46 Cr. If the expansion is slow 
enough that most 45 Cr decays at temperatures 
where 45 V(p, 7) 46 Cr is still rapid, then flow out of 
the N=22 nuclei occurs via (3 + decay out of 46 Cr. 
In this case 46 Ti is synthesized rather than 45 Sc. 

49 Ti originates from the the N=24 nuclide 
49 Mn. At T 9 « 1.4 nuclei with N=24 are divided 
roughly equally between Mn and Fe. Un- 
certainties in the proton separation energies and 
lifetimes of these nuclei are small. 49 Mn does have 
a low lying excited state at 382 keV which is ther- 
mally populated at low temperatures. However, 
49 Mn is a nucleus with Z=N+1 that is expected 
to have ground and excited state decay rates that 
are dominated by super-allowed Fermi transitions 
which are almost independent of excitation energy. 

Lastly, we turn our attention to flow out of the 
N=32 isotones which are progenitors of 60 Zn and 
63 Cu. Proton-rich nucleosynthesis near 64 Ge has 
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been extensively discussed in the X-Ray Burst lit- 
erature (e.g. Brown et al. 2002). Uncertainties in 
basic nuclear properties important for synthesis of 
60 Zn arc small. This is not true for 63 Cu, which is 
formed directly by the decay of 63 Ga. 63 Ga has a 
J n = (5/2)"excited state at 75.4 keV which dom- 
inates the partition function at Tg w 1.5 since the 
ground state has J = 3/2. The weak lifetime of 
this excited state is experimentally undetermined 
(as are the weak lifetimes of all short lived excited 
states) and could easily be a factor of five longer or 
shorter than the quite long ground state lifetime 
of ~ 32 sec. This translates into an uncertainty of 
a factor of several in the inferred 63 Cu yield. 

The influence of possible uncertainties in the 
time scale, entropy, and electron fraction charac- 
terizing the different trajectories can be seen from 
the results in Tabic 2. Modest changes in the out- 
flow parameters result in factors of ~ 2 changes in 
yields of the most important isotopes. This is ev- 
ident by the quite different efficiencies with which 
the lower entropy bubble and higher entropy wind 
synthesize 45 Sc and 49 Ti. 

So far we have not considered the influence of 
neutrino interactions, except implicitly through 
the setting of Y c . If matter remains close to 
the neutron star, neutrino capture and neutrino- 
induced spallation may compete with positron de- 
cay, even on a dynamic time scale. However, neu- 
trino capture alone cannot act to accelerate nu- 
clear flow past waiting point nuclei and allow syn- 
thesis of the heavier proton-rich elements. The 
reason is that the neutrino capture rates on the 
waiting point nuclei arc about the same as the 
rate of neutrino capture on a free proton (Wooslcy 
et al. 1990). Every capture of a neutrino by a 
heavy nucleus is accompanied by a capture onto a 
free proton. The electron fraction is then rapidly 
driven to 1/2 since the neutron produced in this 
way immediately goes into the formation of an a- 
particle. This is analogous to the "a-effect" dis- 
cussed in the context of late-time winds (Fuller & 
Meyer 1995; Meyer et al. 1998). 

4. Conclusions and Implications 

The important news is that, unlike simulations 
of a few years ago, there is no poisonous overpro- 
duction of neutron-rich nuclei in the vicinity of 
the N = 50 closed shell (Woosley et al. 1994). 



When followed in more detail (i.e. mainly with a 
better, spectral treatment of the neutrino trans- 
port), weak interactions in the hot convective bub- 
ble drive Y e back to 0.5 and above so that most of 
the mass comes out as 56 Ni and 4 He. Since 56 Fe 
and helium are abundant in nature, this poses no 
problem. 

Beyond this it is also interesting that the 
proton-rich environment of the hot convective 
bubble and early neutrino-driven wind can syn- 
thesize interesting amounts of some comparatively 
rare intermediate mass elements. If the total mass 
of SN ejecta with Y c > 0.5 is larger than a few 
hundredths of a solar mass, these proton-rich out- 
flows may be responsible for a significant fraction 
of the solar abundances of 45 Sc, 64 Zn, and some 
Ti isotopes, especially 49 Ti. 

However, these ejecta do not appear to be im- 
plicated in the synthesis of elements that do not 
have other known astrophysical production sites. 
For example, Sc can be produced explosively, while 
64 Zn can be synthesized in a slightly neutron-rich 
wind. It seems unlikely that consideration of nu- 
cleosynthesis in proton-rich outflows will lead to 
meaningful constraints on conditions during the 
early SN. 

Since the conditions in the hot convective bub- 
ble resemble in some ways those of Type I X-ray 
bursts (high temperature and proton mass frac- 
tion), we initially hoped that the nuclear flows 
would go higher, perhaps producing the p-process 
isotopes of Mo and Ru. Such species have proven 
difficult to produce elsewhere and the rp-process 
in X-ray bursts can go up as high as tellurium 
(Schatz et al. 2001). Unfortunately the density 
is much less here than in the neutron star and 
the time scale shorter. Proton-induced flows are 
weaker and the leakage through critical waiting 
point nuclei is smaller. Using the present nuclear 
physics, significant production above A = 64 is un- 
likely. However, heavier nuclei can be produced in 
ejecta that are right next to these zones but with 
values of Y c considerably less than 0.50 (Hoffman 
et al. 1996). 

This work was performed under the auspices of 
the U.S. Department of Energy by University of 
California Lawrence Livermore Laboratory under 
contract W-7405-ENG-48. HTJ enjoyed discus- 
sions with Matthias Licbcndorfcr. RB and HTJ 
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uation and visualization, and acknowledge sup- 
port by the Sonderforschungsbereich 375 "Astro- 
Particle Physics" of the Deutsche Forschungsge- 
meinschaft. The supernova simulations were done 
at the Rechenzentrum Garching and at the John 
von Neumann Institute for Computing (NIC) in 
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Table 1 

Reaction Network Used for the Present Calculations 



Element 


N min a 


N b 

^ max 


Element 


N ■ a 

i ' mm 


N b 

i ^ max 


Element 


^ mm 


N b 

iN max 


H 


1 


2 


He 


1 


4 


Li 


3 


6 


Be 


3 


8 


B 


3 


9 


C 


3 


12 


N 


4 


14 


O 


5 


14 


F 


5 


17 


No 


6 


21 


Na 


6 


33 


Mg 


6 


35 


Al 


7 


38 


Si 


8 


40 


P 


8 


42 


S 


8 


44 


CI 


8 


46 


Ar 


9 


49 


K 


11 


51 


Ca 


10 


53 


Sc 


13 


55 


Ti 


12 


58 


V 


15 


60 


Cr 


14 


62 


Mn 


17 


64 


Fc 


16 


66 


Co 


19 


69 


Ni 


18 


71 


Cu 


21 


73 


Zn 


21 


75 


Ga 


24 


77 


Go 


23 


80 


As 


26 


82 


So 


25 


84 


Br 


28 


86 


Kr 


27 


88 


Rb 


31 


91 


Sr 


30 


93 


Y 


33 


95 


Zr 


32 


97 


Nb 


35 


99 


Mo 


35 


102 


Tc 


38 


104 


Ru 


37 


106 


Rh 


40 


108 


Pd 


40 


110 


Ag 


41 


113 


Cd 


42 


115 


In 


43 


117 


Sn 


44 


119 


Sb 


46 


120 


Te 


47 


121 















a Minimum neutron number included for the given clement. 
b Maximum neutron number included for the given element. 
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Table 2 

Characteristics of nucleosynthesis in some representative bubble trajectories 



Trajectory 




s/k b 


Thorn (SCC) 


m/M 


Production (T dyn = lom)* 


Production (Kepler Based Extrapolation) 4 * 


1 


0.500 


18.4 


0.086 


9.25c-04 


59 Co(0.33) 


64 Zn(0.17) 












64 Zn(0.12) 


59 Co(0.17) 












49 Ti(0.10) 


49 Ti(0.16) 


5 


0.502 


15.9 


0.066 


7.05c-04 


59 Co(0.17) 


64 Zn(0.30) 












63 Cu(0.15) 


49 Ti(0.14) 












49 Ti(0.12) 


60 Ni(0.09) 


10 


0.505 


21.7 


0.062 


3.58c-04 


59 Co(0.07) 


64 Zn(0.10) 












63 Cu(0.05) 


49 Ti(0.09) 












49 Ti(0.05) 


46 Ti(0.04) 


20 


0.513 


17.8 


0.104 


4.63c-04 


45 Sc(0.14) 


49 Ti(0.19) 












46 Ti(0.06) 


64 Zn(0.07) 












42 Ca(0.05) 


60 Ni(0.05) 


30 


0.521 


26.2 


0.047 


2.67c-04 


59 Co(0.03) 


49 Ti(0.08) 












45 Sc(0.02) 


64 Zn(0.03) 












63 Cu(0.02) 


60 Ni(0.02) 


35 


0.524 


26.9 


0.062 


2.28c-04 


45 Sc(0.07) 


49 Ti(0.13) 












42 Ca(0.04) 


46 Ti(0.03) 












46 Ti(0.03) 


64 Zn(0.03) 


40 


0.545 


40.6 


0.024 


3.12c-04 


42 Ca(0.04) 


49 Ti(0.25) 












45 Sc(0.04) 


46 Ti(0.06) 












46 Ti(0.03) 


45 Sc(0.04) 



a Listcd here are the three nuclei with the largest production factors. The production factor for each nucleus is given in parenthesis 
next to the nucleus. 
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Table 3 

Production factors in the early wind 



Y e 


s/kb 


T"i 1 ^r>c 1 
' horn 




r iuu ul Liuii ^ ; dyn — ' horn ) 


P roduction ( rvcplcr B ascd Ext rapolation) a 


0.551 


54.8 


0.131 


1.53o-03 


45 Sc(1.73) 


49 Ti(2.02) 










49 Ti(0.97) 


46 Ti(0.70) 










46 Ti(0.87) 


45 Sc(0.36) 


0.558 


58.0 


0.127 


6.40O-04 


45 Sc(0.95) 


49 Ti(1.09) 










49 Ti(0.52) 


46 Ti(0.38) 










46 Ti(0.48) 


45 Sc(0.20) 


0.559 


76.7 


0.099 


6.80o-04 


45 Sc(0.60) 


49 Ti(1.07) 










49 Ti(0.38) 


46 Ti(0.41) 










46 Ti(0.31) 


45 Sc(0.22) 


0.560 


71.0 


0.112 


4.80o-04 


45 Sc(0.55) 


49 Ti(0.79) 










49 Ti(0.31) 


46 Ti(0.29) 










46 Ti(0.27) 


45 Sc(0.15) 


0.568 


74.9 


0.059 


8.00O-04 


45 Sc(0.55) 


49 Ti(1.25) 










46 Ti(0.35) 


46 Ti(0.47) 










49 Ti(0.35) 


45 Sc(0.25) 


0.570 


76.9 


0.034 


1.04o-03 


46 Ti(0.38) 


49 Ti(1.49) 










45 Sc(0.35) 


46 Ti(0.57) 










42 Ca(0.31) 


45 Sc(0.31) 



a Listed here are the three nuclei with the largest production factors. The production factor for each nucleus is given 
in parenthesis next to the nucleus. 



Table 4 

Production factors integrated over the different bubble trajectories 



nucleus 


Production (T dy „ = r ho m) 


Production (Kepler Based Extrapolation) 


Major Progcnitor(s) 


59 Co 


2.81 


0.37 


59 Cu. 59 Zn 


49 T ; 


2.00 


6.53 


49 Mn 


63 Cu 


1.91 


0.28 


63 Ga, 63 Gc 


45 Sc 


1.65 


1.33 


45 Cr 


64 Zn 


1.28 


3.61 


64 Gc 


46 T ; 


1.22 


1.97 


46 Cr 


60 Ni 


1.10 


1.81 


60 Zn 


42 Ca 


1.04 


0.46 


42 Ti 



Table 5 

Integrated production factors for the early wind 



nucleus 


Production (T dy „ = T ho m) 


Production (Kepler Based Extrapolation) 


46 Sc 


4.74 


1.50 


49 Ti 


2.83 


7.70 


46 Ti 


2.66 


2.81 


42 Ca 


2.16 


0.46 


Sly 


1.09 


0.90 


50 Cr 


0.56 


0.09 
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Table 6 

Characteristics of Nucleosynthesis in Neutron Rich Trajectories 





m/M Q 


Production" 


0.470 


6.40o-05 


74 Sc(6.59) 






78 Kr(4.25) 






64 Zn(1.35) 


0.475 


7.98o-05 


64 Zn(1.36) 






74 Se(0.85) 






78 Kr(0.78) 


0.480 


1.59o-04 


64 Zn(1.49) 






78 Kr(0.34) 






68 Zn(0.30) 


0.485 


3.36o-04 


62 Ni(0.92) 






58 Ni(0.35) 






64 Zn(0.23) 


0.490 


6.24o-04 


62 Ni(1.21) 






58 Ni(0.42) 






66 Zn(0.13) 


0.495 


1.36o-03 


62 Ni(1.30) 






58 Ni(0.41) 






61 Ni(0.23) 



a Listed here are the three nuclei 
with the largest production factors. 
The production factor for each nu- 
cleus is given in parenthesis next to 
the nucleus. 
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Table 7 

Nuclei for which thermal weak rates are included 



Atomic mass Elements 11 



21 


F, Mg, Na, Ne, O 


22 


Mg, Na, Nc 


23 


F, Mg, Na, Nc 


24 


Mg, Na, Nc, Si 


25 


Mg, Na, Ne, Si 


26 


Mg, Na, Si 


27 


Mg, Na, P, Si 


28 


Mg, Na, P, S, Si 


29 


Mg, INa, r, D, ol 


30 


P, S, Si 


31 


CI, P, S, Si 


32 


CI, P, S, Si 


33 


CI, P, S, Si 


34 


CI, P, S, Si 


35 


CI, K, P, S 


36 


Ca, CI, K, S 


37 


Ca, CI, K, S 


38 


Ca, CI, K, S 


39 


Ca, CI, K 


40 


Ca, CI, K, Sc, Ti 


41 


Ca, CI, K, Sc, Ti 


42 


Ca, K, Sc, Ti 


43 


Ca, CI, K, Sc, Ti 


44 


Ca, K, Sc, Ti, V 


45 


Cr, K, Sc, Ti, V 


46 


Cr, K, Sc, Ti, V 


47 


Cr, K, Sc, Ti, V 


48 


Cr, K, Sc, Ti, V 


49 


Cr, Fe, K, Mn, Sc, Ti, V 


50 


Cr, Mn, Sc, Ti, V 


51 


Mn, Sc, Ti, V 


52 


Fe, Mn, Ti, V 


53 


Cr, Fe, Mn, Ti, V 


54 


Cr, Fe, Mn, V 


55 


Cr, Fc, Mn, Ti, V 


56 


Cr, Fc, Mn, Ni, Sc, Ti, V 


57 


Cr, Cu, Fe, Mn, Ni, Ti, V, Zn 


58 


Cr, Cu, Fe, Mn, Ni, Ti, V 


59 


Cr, Cu, Fc, Mn, Ni, V 


60 


Cr, Cu, Fe, Mn, Ni, Ti, V, Zn 



a All elements of the given mass for which the 
Fuller, Fowler, Sz Newman (1982) rates were in- 
cluded. 
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200 400 600 

t pb [ms] 

Fig. 1. — Radius versus post-bounce time for se- 
lected "mass shells" in the exploding SN model 
of Janka, Buras, & Rampp (2003), spaced with 
intervals between 0.001 M Q and 0.1 M Q . The sim- 
ulation was carried out in two dimensions until 
about 470 ms after bounce, and was continued in 
spherical symmetry thereafter. During the 2D 
period the lines do not trace the trajectories of 
Lagrangian mass elements but correspond to the 
radii of spheres enclosing certain values of the in- 
tegrated rest mass. The expanding dashed lines 
mark the positions of forward shock and wind ter- 
mination shock, respectively, and the thin dashed, 
wiggled line represents the angularly averaged po- 
sition of the gain radius. The neutrinosphere po- 
sitions of electron neutrinos (solid), electron an- 
tineutrinos (dashed), and muon and tau neutri- 
nos and antineutrinos (dash-dotted) are also indi- 
cated. The explosion sets in when the shock passes 
the infalling interface between the Si layer and the 
oxygen-enriched Si layer (given by the trajectory 
marked with dots) at which the density begins to 
drop steeply and the entropy increases from about 
2.5 to nearly 5 kb per nuclcon. 
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Fig. 2. — Ejecta mass versus Y e of neutrino-heated 
and -processed matter during the convective phase 
until ^470 ms post bounce. The insert shows the 
region around Y c ~ 0.5 in higher resolution. The 
grey shading indicates estimated errors due to the 
limited spatial resolution of the two-dimensional 
simulation (for details, see Buras et al. 2004). 
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Fig. 3. — Density, temperature, Y c , and entropy as functions of post-bounce time along the trajectories of 
mass elements around an enclosed baryonic mass of 1.41 M©. The elements first follow the rise of temperature 
and density in the outer layers of the contracting neutron star and then enter a phase of very rapid expansion 
when they are ejected in the neutrino-driven wind. The curves are labeled by the time the mass elements 
cross a radius of 100 km. The collision with the slower preceding ejecta occurs through a wind termination 
shock and is visible as a non-monotonicity of the density and temperature, associated with an entropy 
increase of 10-15 kb per nucleon. 
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Fig. 4. — Convective bubble around the central 
neutron star in the two-dimensional SN model of 
Janka, Buras, & Rampp (2003). Time in this fig- 
ure is measured from the onset of collapse, core 
bounce was at 175 ms, the length scale is in km. 
The panels on the left show the temperature in 
Kelvin, those on the right Y c . The upper plots 
give radial profiles for all angular bins of the po- 
lar coordinate grid of the simulation, which was 
carried out in a lateral wedge of ±43.2° (with pe- 
riodic boundary conditions) around the equatorial 
plane. The latter is indicated by white diagonal 
lines. The positions of the tracer particles at the 
onset of the explosion are marked by crosses in 
the lower panels. Their positions were chosen (by 
post-processing the finished simulation) such that 
the Y e distribution of the final ejecta was appro- 
priately represented. 



20 




Fig. 5. — Illustration of the density evolution for 
hot-bubble trajectory 10 in the two-dimensional 
SN model of Janka, Buras, & Rampp (2003) (left 
panel) and as extrapolated for the present nucle- 
osynthesis calculations (right panel). 
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Fig. 6. — Evolution of temperature in the trajec- 
tory displayed in Fig. 5. 
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Fig. 7. — Overproduction plot for nucleosynthe- 
sis in some of the tracer particle trajectories in 
the convective bubble of the SN model of Janka, 
Buras, & Rampp (2003). Results for two assump- 
tions about the dynamic time scale are shown for 
each trajectory. 
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Fig. 8. — Final proton and a particle mass frac- 
tions as a function of Y c for the different Kepler 
extrapolated trajectories. The line labeled with 
X p (theory) is the proton mass fraction calculated 
under the assumption that all available nucleons 
are locked into alpha particles or other nuclei with 
equal and even numbers of protons and neutrons. 
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Fig. 9. — Evolution of nuclear abundances as a 
function of neutron number in the Kepler based 
extrapolation for the hot-bubble trajectory 26 
(Y e = 0.5172). 
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